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Abstract 

We consider the drift of a polymer chain in a disordered medium, which 
is caused by a constant force applied to the one end of the polymer, under 
neglecting the thermal fluctuations. In the lowest order of the perturbation 
theory we have computed the transversal fluctuations of the centre of mass 
of the polymer, the transversal and the longitudinal size of the polymer, and 
the average velocity of the polymer. The corrections to the quantities under 
consideration, which are due to the interplay between the motion and the 
quenched forces, are controlled by the driving force and the degree of poly- 
merization. The transversal fluctuations of the Brownian particle and of the 
centre of mass of the polymer are obtained to be diffusive. The transversal 
fluctuations studied in the present Letter may also be of relevance for the 
related problem of the drift of a directed polymer in disordered media and its 
applications. 

PACS numbers: 05.40.+j., 36.20C, 61.41. +e. 
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The behaviour of a polymer chain in disordered media is of interest for different applica- 
tions such as gel electrophoresis of biopolymers (see [[I]- || and references therein). It is also 
expected that the study of the behavior of a polymer in a disordered medium may be useful 
for a better understanding such a fundamental problem as reptation 0- ||. Despite the large 
interest in this problem in the last decade 0- ]T3| no complete quantitative understanding 
of the behaviour of the polymer chain has been achieved so far. The reason is the quenched 
nature of the disorder, which after the average results in an effective attraction between the 
polymer segments and is responsible that the problem under consideration becomes a strong 
coupling problem. The standard theoretical tools such as perturbation theories and also the 
renormalization group method break down in a number of interesting cases for these kind 
of problems due to the runaway of the effective coupling constant. Thus, the polymer in a 
disordered media is characterized by a nontrivial interplay between the thermal effects and 
the disorder, which so far is beyond of a quantitative analytical treatment. The interplay 
between disorder and thermal effects is also of great importance for different kind of creep 
phenomena [IB|- |I7| (and references therein). 



In the present Letter we will present the results of the study of the simpler problem of 
the drift of a polymer chain in disordered media under the action of a constant force applied 
to one end of the polymer by neglecting the thermal effects. The motivation of the present 
problem is the electrophoresis of biopolymer molecules in porous media. The drift of a 
polymer in disordered media can be considered as a further example of the driven dynamics 
of manifolds (flux lines, interfaces, etc.) in disordered media, where a considerable progress 



was achieved in last years ||18|| . One of the interesting aspects of the driven dynamics of 
the polymer chain in a disordered medium is the generation of the transversal fluctuations, 
which is due to the interplay between the motion and the quenched random forces. This 
effect, which to our knowledge have not been studied previously, appears already in the 
case of a Brownian particle (at zero temperature). We will argue in this Letter that the 
transversal fluctuations of the centre of mass of the polymer chain (and of the Brownian 
particle) generated due to the interplay between the disorder and motion are diffusive. We 



expect that such transversal fluctuations may be of interest for directed manifolds (flux lines, 
crack propagation, etc.). 

We start the consideration with the Langevin equation for a discrete polymer chain || 

Ot ■ 

f-j£ = M r i+i + r i-i ~ 2r i) + F i( r i) + Vh (!) 

where rf give the positions of monomers (i = 1, N) in ci-dimensional space (a = 1, d), 
f is the monomer friction coefficient. The elastic forces are associated with the potential 
U — |r Sili( r « ~ r i-i) 2 - Notice that the elasticity constant x is given by x = dkT/l 2 , I is 
the statistical segment length, and iV = L/l is the number of monomers. In the following 
we will consider the drift of the polymer chain at zero temperature, so that we will consider 
k as a constant, which is temperature independent. The force -Fj(r) = F5n + g(r) with F 
being the constant driving force acting on the 1st segment of the polymer, and g(r) is the 
random quenched force. We assume that the latter is Gaussian distributed with the zero 
mean and the correlator 

[<f(r)/(r')] =C aP {r-r'). (2) 

The Fourier transform of C al3 (r), C a @(q), is assumed to decay exponentially for large q with 
a width a. In particular C aa (q) has the shape C aa (q) = Aexp(— a 2 q 2 ) with a = x,y,z in 
3d. The thermal forces are Gaussian distributed with the correlator given by 

< rj?(t) V ?(t') > = 25 a %fkT5(t - t'). (3) 

The Langevin equation (|]) can be formulated in a standard way by using the path inte- 
grals, which in the limit of a continuous chain || results in the following expression for the 
transition probability density for a polymer configuration r(s,t) 

P(r(s),t; r°(s), t Q ) = [ ( ^ () Dr(s 7 1) [ Dp( Sj t) exp(-S), (4) 

Jr(s,t )=r°(s) J 

where the action for the continuous polymer chain reads 
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S = D/i 2 I ds f dt'p(s,t') 2 -i I ds f dt' P {s,t'){ dr ^P -7/iVHM') ~ Ms{r{s,t')) 

JO Jt Jo Jt 

(5) 



with D = fkT/l, /i 1 = f/l, 7 = xZ, and / being the statistical seg- 
ment length. The generating functional of the non interacting polymer, Zo(j,j) = 
L(s,t )=r°(s) Dr ( s > *) / D P( S > f ) ex P(-^ + i J L ds /* dt'(r(s, t')j(s, f ) +p(s, t')j(s, t')), is given 
by 

Z (j J) = exp(-i / dh I dhjk^GHh - t 2 )~j k {t 2 ) + % j dt l3k (t 1 )Gl(t 1 -t 2 )Zl + 

J to J to Jto 

rt 

dti I dt 2 jk(t 1 )G° k (ti - t 2 )ng k - 

Jto 



to J to 

2 r t 



dti / ^2j fc (ti)(G^(ti-t 2 )+G^(t 2 -t 1 )-G^(t 1 +t 2 -2t ))j fc (t 2 ), (6) 

to J to 



+ 



where the Fourier transform of a quantity a(s) is defined as follows a(s) = YlT=o QskCtk with 
Q sk = y/2/L cos(tt sk/L), and Q s0 = 1/VL, G° k (t) = exp(-~ffik 2 t)9(t) with 0(t) being the 
Heaviside function, and £ k is the Fourier transform of r(s,t = 0). The last term in the 
exponential of Eq.(||) disappears for zero temperatures. The generating functional Z(j,j = 
0) at zero temperature can be written in the symbolic form as 

Z(j) = e W (f dh f dt 2 I ' ds, ( L ds 2 j3t{ti)Gl{t-t l )Q Slk C^{q)Q S2k Gl{t-t 2 ) ] l{t 2 ^ 

JO JO Jo Jo Jq 

if dh [ dt 2 j k (t 1 )G° k (t 1 -t 2 )ixg k ), (7) 

Jt J t 

where the source j k {t) in Eq.(|7|) has to be modified in each order of the perturbation ex- 
pansion as follows: j k (t) -> j k (t) + qQ Slk S(t - tt) - qQ S2k S(t - t 2 ), where q, s x , s 2 , h, 
and t 2 are variables associated with the disorder correlator, in the 1st order, jk{t) — > 
Jk(t) + qiQ Slk S(t - ti) - qiQ S2k b{t - t 2 ) + q 2 Q S3k 8(t - t 3 ) - q 2 Q Sik 6(t - £ 4 ) in the 2nd order 
and so on. 

We now will study the behaviour of the ideal polymer chain under the action of the 
force acting at one end of the polymer under neglecting thermal fluctuations. Due to the 
modification of jk(t) the last term in (|7|) gives rise (in the steady state limit) to the sum 



X] (Qaik ~ Q s 2 k) / {itk / L) 2 = y/L/2(s2 — Si)(l — (si + s 2 )/2L). Both the transversal fluctua- 

k=l 

tions of the centre of mass and the transversal size of the polymer can be expressed by the 
quantity 

= ^Z(j)/5jl(t)5f k2 (t). (8) 

The transversal size of the polymer is obtained from Eq.(|8|) with k\ and hi nonzero and a, 
(3 7^ z. The analytical expression associated with Eq.(^) reads 

fi 2 r 2 [ dsi [ ds 2 [ dt 2 ^ dU f C af '(q)Q Slkl Q S2k2 Gl 2 (t-t 2 )G ki (t-t 1 ) 

JO JO Jo Jo Jq 

exp(ig z 5( fc (ti - t 2 ) + iq'{s 2 - s x )(l - (s t + s 2 )/2L)), (9) 

where g k = fJ-F/L and q' = q z F/ (V27). Due to the fact that we are interested in large t, we 
have neglected the time dependent part associated with the second term in the exponential 
of Eq. @ . Integrations over times in @ can be performed straightforwardly. The integrals 
over si and s 2 in (|j) can be also computed analytically. It appears that the main contribution 
in the latter originates from k\ = k 2 (the term with k\ 7^ k 2 vanishes for large F), so that 
we obtain 

A(k,q)= / dsi / ds 2 cos(fcsi) cos(fcs 2 ) exp(iq(s 2 — Si)(l — (si + s 2 )/2L)) = 
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(nL/Aq)((C(^/L^Tq(q + k)) + C{^hJVq{q - k)) 2 + 

{S(y/Lfa{q + A;)) + ^(V 7 ^^ - A:)) 2 ), (10) 

where C(x) and S(x) are Fresnel functions. In the limit of large L the function A(k,q) 
tends to the function (irL/2q)(6(\ q \ — | fc |). After performing the integrations over ^ 
and Sj (i = 1,2) in Eq.(|8|) and using the property of A(k,q) for large L, the mean-square 
transversal size of the polymer chain is obtained to the first order in disorder strength as 

[Ary = {d -^f L \ ^A/a d + 0(A 2 ). (11) 

The fluctuations of the centre of mass of the polymer chain to the first order in disorder 
strength can be studied by using (§|f)|) with k\ = k 2 = 0. The mean-square transversal 



displacement of the centre of mass of the polymer chain is obtained after performing the 
integrations over times ti, t 2 and si, s 2 as 

[ArKt)} = {d ~^ tL ^ d -'A/a d + 0(A 2 ) (12) 

Setting L = I in Eq.fJT^) yields the mean-square transversal displacement for a Brownian 
particle. 

We now will derive Eq. fll2|) for a Brownian particle in a qualitative way. Due to the 
motion in the direction of the force F, the particle experiences the random force in the 
transversal directions too. The transversal fluctuations can be estimated according to the 
Langevin equation of the particle as f 2 x 2 r /t 2 ~ 5F 2 with 5F 2 ~ aA/a d /lp being the square 
of the typical value of the random force and lp = fiFt being the characteristic length. 
It follows from Eq.(|T2|) and from the above qualitative consideration that the transversal 
behaviour of the centre of mass of the polymer chain and of the Brownian particle is diffusive. 
The analysis of higher-order corrections to (|12|) shows that they are linear in t, so that the 
transversal fluctuations of the particle are expected to be diffusive. Rewriting flT2"D in terms of 
the unperturbed velocity v c $ = \iFj L yields that the mean-square transversal displacements 
for both the Brownian particle and the polymer chain coincide. 

Excluding the time t in fll2|) in favor of the arc length L by using the relation fi'jt ~ L 2 , 
which is obtained from the unperturbed response function exp(— / ~ffi(irk/L) 2 t) or from the 
relation v c t ~ 5r z with v c and 5r z given by Eqs. (p!6| , p!8|) P| up to order 0(A°), we find that 



Eq.(|12"D is compatible with Eq.flll]). To our knowledge the transversal fluctuations given by 
( |TT| - |T^) have not been studied before. Notice that the temperature is zero and the transversal 
fluctuations are due to the interplay between the motion and the frozen in random forces. 
To study the average velocity of the polymer we will consider the quantity 

{ti(t)] = r l 8Z{j)/5f k {t) (13) 



The latter argument was suggested by one of the referee. 
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in the limit k = 0. The analytical expression associated with (|13| ) is given by 

[%(t)] = VLg k t + f iH- 2 f2 f {-iq a )C za {q) [ d Sl [ ds 2 [ dt 2 I ' * dt x Q sa0 Q slk Q S2k 

k=1 Jg JO JO Jo Jo 

G° k (t 2 - <i)exp(ig z (/ fc (*i - h) + iq\s 2 - si)(l - (si + s 2 )/2L)). (14) 

After performing in ([L4]) the integrations over times and arc lengths we arrive at 

, r3/2 rn xx (n\ r\/L\<iz\»F/{iivL) 

v c (t) = ^F/L(l - v 7 ^^ / / dx ^ x " + X ) + -)■ ( 15 ) 



9 Vl^l "^0 



The upper limit in the integration over x can be put to infinity for large driving force F, so 
that ( |i~5|) simplifies to 

« c (0 = imF/LQ. - c ^< 2F ^ (A/a d ) + ...), (16) 

with c v = vr 2 /(V2r(3/4)). Eq.(|16D gives the first-order correction to the velocity, so that 
the correction term in ( |16|) is expected to be small. This condition is fulfilled for large 
driving force F. For finite F and large L the disorder drastically affects the behaviour of 
the polymer chain. According to Eq. fll6|) , the validity of the perturbation expansion for a 
given force F is restricted to the arc length of the polymer L, L < F 5 ^ 3 . 

Let us compare Eq.fll^) with the driven motion of an interface in disordered media 
|i~8|| . The expansion parameter of the perturbation expansion (|i~6D, v /toX 3 / 2 / (Z 2 F 5 / 2 )(A/a d ), 
depends on F as F~ 5 ^ 2 instead of F~ 3 ^ 2 as it is the case for interfaces in d = 1. The additional 
factor 1/F for the polymer problem is a consequence of nontrivial integrations over s\ and 
s 2 , which are due to the absence of the factor 5^(xi — x 2 ) in the disorder correlator (^) 
and to the fact that the driving force acts on the extremity of the polymer. The additional 
factor 1/F has the consequence that the expansion parameter cannot be expressed in terms 
of the force density F/N. The latter expresses the nontrivial carrying over of the driving 
force F to the monomers of the polymer chain. 

To study the longitudinal size of the polymer we consider Eq . fll~3|) with k ^ 0. The 
analytical expression associated with the latter is 
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[£(*)] =f% + /i 2 r 2 V f(-tq a )C za (q) [ d Sl [ ds 2 [ dt 2 I' dt x Q S2P Q slk Q S2k 

Kl^V k=1 Jq JO JO JO JO 

G°(t - £i)G°(t 2 - h) exp(iq t g k (h - h) + iq'(s 2 - si)(l - (si + s 2 )/2L)). (17) 



The average value [r 2 (0) — r z (L)] = [Sr z (L)] is obtained in terms of the normal coordinates 



as 2a/2/L Y1T=i tefc-i] • Performing the integrations over the times and arc lengths results 
in 

with c r = 7r 2 /(2r(3/4)). Notice that the disorder weakens the stretching of the polymer. 
This can be understood as follows. It is evident that the effect of the disorder results in a 
decrease of the unperturbed velocity \iFjL. The decrease of the velocity due to disorder 
can be interpreted as the unperturbed motion due to a smaller force, which in virtue is 
accompanied by a less stretching of the polymer. Thus, this argument supports the non 
evident result of Eq.flTBD that the stretching of the polymer in the longitudinal direction will 
be diminished due to disorder. 

To conclude, we have considered the drift of a polymer chain in disordered media under 
the action of a constant force applied to one of the extremities of the polymer by neglecting 
the thermal effects. The dynamics of the polymer is governed by the interplay between the 
motion and the quenched random forces. One of the significant effects of this interplay is 
the generation of transversal fluctuations of the polymer chain. The latter take also place 
for a Brownian particle. We have also found that disorder results in a weakening of the 
elongation of the polymer in the longitudinal direction, which is caused by the friction of 
the polymer with the surrounding medium. Due to the specific dependence of the quantities 
under investigation on the driving force F and the length of the polymer L, the fluctuations 
induced by the drift of the polymer may compete with thermal fluctuations. An interesting 
question is the experimental manifestation of these fluctuations. However, the first-order 
corrections studied in the present Letter do not allow us to make predictions on the total 
effect of the disorder on the quantities under consideration. The only conclusion we may 
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drawn on this matter is that for finite driving force the effect of disorder increases with the 
polymer length. The transversal fluctuations studied in the present Letter may be also of 
relevance for the related problems of the drift of a directed polymer in disordered media, 
crack propagation in disordered media etc. To our knowledge, the transversal fluctuations 
in these problems have not been studied so far. 
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